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I.  INTRODUCTION 


This  final  technical  report  on  Contract  No.  N00014-74-C-0154 
summarizes  a research  Investigation  conducted  by  Desmatlcs,  Inc.  under 
sponsorship  of  the  Office  of  Naval  Research.  The  Desmatlcs  study  has 
been  devoted  to  statistical  research  In  a number  of  problem  areas  associ- 
ated with  the  Impact  acceleration  and  ship  motion  programs  being  conducted 
by  the  Naval  Aerospace  Medical  Research  Laboratory  (NAMRL)  Detachment. 

The  following  section  of  this  report  briefly  reviews  research  results 
In  the  main  problem  areas  addressed  by  Desmatlcs.  The  final  section  pro- 
vides a reference  list  of  technical  reports  and  other  publications  prepared 
under  this  contract. 


II.  RESEARCH  RESULTS 

During  the  technical  effort  performed  under  Contract  No.  N00014-7A-C-0154, 
Desmatlcs  has  conducted  statistical  research  associated  with  NAMRL  Impact 
acceleration  and  ship  motion  programs.  The  following  paragraphs,  which  review 
research  results,  briefly  describe  the  three  primary  problem  areas  considered 
by  Desmatlcs. 

A.  HEAD/NECK  DYNAMIC  RESPONSE  MODEL 

A portion  of  the  Desmatlcs  technical  effort  has  been  devoted  to  develop- 
ment of  a predictive  model  of  dynamic  response  of  the  human  head/neck  system 
to  -G^  Impact  acceleration.  This  model,  based  on  empirical  data  only,  I 

requires  no  assumptions  about  the  mechanical  structure  of  the  system.  It 
predicts  the  time-varying  Inertial  response  of  the  unrestrained  human  head 
to  -G^  Impact  acceleration.  Using  the  linear  acceleration  time  profile  of 
neck  (T^  vertebra)  motion  as  Input  Information,  the  model  produces  time  pro- 
files for  the  accelerations  at  the  head  center  of  gravity. 

Model  development  (l.e.,  estimation  of  all  required  model  parameters) 
was  based  upon  data  from  NAMRL  Impact  acceleration  experiments.  In  these 
experiments,  volunteer  subjects  were  strapped  securely  to  a steel  chair 
rigidly  mounted  on  a sled.  During  sled  acceleration,  the  subject's  motion 
was  measured  by  Instrumentation  which  Included  two  accelerometer  arrays: 
one  located  at  the  mouth,  and  the  other  on  the  posterior  spinous  process 
I of  the  first  thoracic  vertebra  (T^^). 

Model  parameters  were  estimated  using  the  data  from  six  Instrumented 
human  volunteers,  who  were  each  exposed  to  a number  of  -G^  Impact  accelerations. 
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Data  from  a total  of  24  tests  was  used  to  estimate  the  required  parameters. 

The  basic  model,  representing  a group  of  six  Individuals,  was  developed 
using  the  linear  X acceleration  measured  at  vertebra  as  the  Input,  the 
linear  X and  Z accelerations  and  angular  Y acceleration  of  the  head  CG  as 
outputs.  Then,  for  a second  Independent  set  of  24  Impact  runs  the  model  was 
used  to  predict  the  head  accelerations  produced  by  each  T^  acceleration.  In 
each  Instance,  the  model's  response  was  visually  compared  with  the  actual 
response  measured  from  the  subject.  In  general,  the  model  and  measured 
responses  were  In  good  agreement,  approximately  as  good  as  the  agreement 
between  the  responses  of  the  same  subject  undergoing  two  Identical  Impact 
accelerations.  This  "eyeball"  evaluation  was  augmented  with  a statistical 
analysis  of  model  performance  based  on  F-tests  designed  to  test  model  bias. 

From  evaluation  of  model  performance,  which  was  based  on  data  Independent 
of  that  used  In  model  development.  It  was  found  that  variations  between  model 
output  and  measured  head  response  were  comparable  to  variations  between  the 
head  responses  of  subjects  undergoing  the  same  Impact  accelerations.  This 
Indicates  that  the  model  will  be  useful  In  subsequent  studies  of  human  response 

to  Impact  acceleration.  > 

i 

I 

A study  of  model  sensitivity  to  design  changes  suggests  that,  based  on  ) 

j 

the  available  evidence,  the  selected  model  Is  a satisfactory  compromise  | 

! 

between  performance  and  simplicity.  In  addition  to  Its  relatively  good  per- 
formance, the  model  has  the  advantage  that  Its  development  Is  based  only  on 
observed  Inertial  response  data.  Model  accuracy  Is,  therefore,  not  dependent 
on  often  hypothetical  parameters  defining  the  mechanical  structure  of  the 
head/neck  dynamic  system. 

This  aspect  of  the  Desmatlcs  research  has  been  discussed  In  technical 


J 

i 


-3- 


e: 


presentations  [2,5]  at  two  scientific  meetings  and  in  a journal  article  [8]. 

B.  IMPACT  ACCELERATION  INJURY  PREDICTION  MODEL 

The  Desmatlcs  research  effort  has  also  be^n  devoted  to  investigation 
of  head/neck  injury  prediction  models  for  the  impact  acceleration  environment. 
The  general  framework  of  a statistical  model  has  been  discussed  in  two  tech- 
nical reports  [1,3].  The  resulting  model  may  be  used  to  predict  the  probability 
of  head/neck  Injury.  As  a basis  for  discussing  the  model,  consider  the  Impact 
acceleration  situation  in  which  the  torso  is  well-restrained,  but  the  head 
and  neck  are  unrestrained.  In  this  situation  the  problem  is  one  of  predicting 
whether  a human  of  given  anthropometric  characteristics  will  sustain  Injury 
if  exposed  to  Impact  acceleration  which  results  in  given  dynamic  response  of 
the  head/neck  system.  A number  of  difficulties  must  be  overcome  to  develop 
an  injury  prediction  model  from  empirical  data.  If  enough  Instrumented  human 
subjects  were  available,  and  could  be  subjected  to  various  acceleration  time 
traces,  a reasonable  prediction  model  would  eventually  result.  Of  course, 
this  procedure  is  not  possible — human  subjects  cannot  be  purposely  Injured. 

In  any  event,  it  must  be  realized  that  the  situation  is  not  deterministic. 
For  example,  even  with  a restrained  torso,  the  same  impact  acceleration  will 
result  in  different  head/neck  response  (for  example,  because  of  initial  head 
position),  and  even  apparently  identical  head  response  for  the  same  person 
may  result  in  Injury  sometimes  and  not  at  other  times.  This  binomial  aspect 
of  Injury  occurrence  defines  a discrete  random  variable  which  must  be  considered. 
To  further  complicate  matters,  the  acceleration  and  dynamic  response  data 
under  consideration  is  time  trace  data. 

Because  dealing  with  the  complete  acceleration  time  traces  of  the  head 


Is  an  Impossible  analytic  task,  a set  of  univariate  head  dynamic  response 
variables  which  may  be  expected  to  be  related  to  Injury  should  be  considered. 
Likely  candidates  Include,  for  example,  linear  and  angular  velocities  and 
accelerations  (average  or  peak) . Although  a number  of  anthropometric  vari- 
ables can  also  be  postulated.  It  Is  probably  reasonable  to  assume  that  their 
effect  within  a species  will  be  minor  when  compared  to  that  of  the  dynamic 
response  variables.  Thus,  It  Is  suggested  that  only  these  latter  variables 
be  considered  In  Initial  model  development.  At  a later  stage,  the  anthropo- 
metric variables  should  prove  Important  In  scaling. 

In  general,  then,  a set  of  k variables,  which  will  be  denoted  by 
X =»  (x^^, . . . ,Xj^) , is  being  considered.  It  is  postulated  that  the  probability 
of  injury  is  some  (unknown)  function  of  these  variables.  Furthermore,  although 
the  function  Is  unknown.  It  will  be  near  zero  In  part  of  x~space,  near  one 
In  another  part,  and  will  Increase  from  near  zero  to  near  one  over  an  inter- 
mediate part.  Experimentally  what  Is  observed  In  a given  situation  Is  only 
an  estimated  value  (either  0 or  1)  of  the  true  probability. 

In  summary,  the  probability  of  Injury  Is  being  considered  as  a function 
of  X.  Thus,  this  probability  may  be  denoted  by: 


P » P (x)  » P (Xj^,  . . ..Xj^) 


Furthermore,  the  observed  value  of  F will  be  denoted  by  y,  where: 


1 If  an  Injury  Is  sustained 
0 If  no  Injury  Is  sustained. 


It  will  be  assumed  that  a logistic  function  provides  a reasonable 
approximation  to  the  function  defining  probability.  The  logistic  function 
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is  given  by: 

^ -1 

P(x)  = {l  + exp[-(gQ  + Eg^x^)]} 

When  this  function  is  used,  all  predicted  probabilities  are  restricted  to 
the  range  (0,1).  Furthermore,  this  function  satisfies  the  conditions  of 
being  near  zero  in  a part  of  x-space,  near  one  in  another  part,  and  increas- 
ing from  near  zero  to  near  one  over  an  intermediate  part.  It  is  also  tractable 
computationally . 

From  a set  of  observed  data,  the  coefficients  (i.e.,  g^,...,gj^)  of 

the  logistic  model  may  be  estimated.  The  estimation  process  is  fairly  complex, 
involving  an  iterative  procedure  which  provides  the  maximum  likelihood  estimates 
This  does  not  pose  an  Insurmountable  problem,  however,  since  the  computer  is 
available. 

Nonetheless,  the  data  input  to  a model  of  this  kind  is  of  necessity 
dichotomous,  requiring  the  use  of  larger  samples  than  required  to  obtain  a 
desired  degree  of  predictive  accuracy  if  the  data  were  continuous.  Thus,  it 
is  of  central  concern  to  investigate  the  degree  of  accuracy  which  may  be 
expected  for  predictions  derived  from  the  model,  and  to  examine  the  sensitivity 
of  such  predictions  to  sample  size. 

A Monte  Carlo  simulation  study  was  undertaken  to  provide  information 
relating  accuracy  to  sample  size  for  selected  model  configurations.  Two 
specific  sets  of  model  parameters  were  considered,  Monte  Carlo  samples  of 
various  sizes  were  generated  for  each,  and  the  accuracy  of  the  resultant  pre- 
dictions were  evaluated  with  respect  to  the  true  probabilities.  This  study 
is  described  in  a Desmatlcs  technical  report  [9]. 

Research  has  also  been  focused  on  applications  to  real  data  and  on 
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statistical  Inferences  based  on  model  results.  A recent  technical  report  j 

[10]  and  technical  presentation  [7]  considered  application  of  a logistic 

model  to  observed  data  from  a set  of  twenty-eight  -G^  accelerator  runs  4 

Involving  subhuman  primates  (Rhesus  monkeys)  with  securely  restrained  torso 

and  unrestrained  head.  The  data  was  collected  by  the  Naval  Aerospace 

Medical  Research  Laboratory  (MAMRL)  Detachment  as  part  of  Its  research  I 

effort  on  acceleration  Impact  Injury  prevention. 

Two  prediction  models  were  constructed  from  this  data,  one  based  on 
head  dynamic  response  only  and  the  other  based  on  sled  acceleration  profile 
only.  The  first  model  was  derived  from  three  variables  distilled  from 
head  dynamic  response  time  trace  data.  These  three  variables  were; 


(1)  peak  head  angular  acceleration  (resultant)  measured  In 
radians /sec 2, 

(2)  peak  head  linear  acceleration  (resultant)  measured  In 
meters/sec2, 

and  (3)  peak  head  angular  velocity  (resultant)  measured  In  ! 

radlans/sec. 

The  second  model  was  based  on  two  variables  describing  sled  acceleration: 

(1)  peak  sled  acceleration  measured  In  G's 

t 

and  (2)  rate  of  sled  acceleration  onset  measured  In  G/sec. 

i 

Surprisingly,  the  model  based  on  sled  profile  variables  provided  || 

better  predictions  than  the  model  based  on  head  dynamic  response  variables.  I 

Because  of  this,  both  sets  of  variables  were  combined  Into  one  overall 
five-variable  set,  which  was  then  used  In  development  of  a prediction  model. 

It  was  found  that,  for  the  relatively  small  amount  of  data  available. 

Inclusion  of  peak  sled  acceleration  and  any  one  of  the  other  variables  ! 
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resulted  in  a model  which  yielded  predictions  in  almost  perfect  agreement 
with  the  observations. 

Since  it  is  intuitive  that  a model  based  on  head  dynamic  response 
should  provide  predictions  which  are  as  good  as  or  better  than  those  from 
a model  based  on  nled  profile,  some  explanation  is  required.  There  are  a 
few  possible  reasons  for  this  anomalous  result.  It  may  be  that  the  wrong 
variables  were  extracted  from  the  head  dynamic  response  time  traces,  and 
other  variables  would  have  more  correctly  conveyed  the  information  within 
these  time  traces.  On  the  other  hand,  the  correct  variables  may  have  been 
selected  but  errors  may  have  been  present  in  their  measurement.  In  addition. 


it  is  possible  that  the  small  sample  size  resulted  in  a spurious  result. 

It  must  be  realized,  however,  that  consideration  of  any  of  the  three 
head  dynamic  response  variables  in  conjunction  with  peak  sled  acceleration 
provides  a perfectly  fitting  model.  Thus,  there  is  evidence  that  head 

|i 

dynamic  response  variables  and  sled  profile  variables  may  be  used  together 
to  provide  good  results.  Nonetheless,  it  is  still  an  open  question  of  why 
head  dynamic  response  variables  alone  did  not  provide  a model  which  performed 
as  well  as  the  one  based  on  sled  profile  variables  alone.  Further  research 
on  this  topic  shou.J  prove  valuable. 

Additional  research  has  been  concerned  with  statistical  inferences  about 
the  parameter  vector  of  the  logistic  model,  and  P(x),  the  true  probability 


of  injury  corresponding  to  acceleration  variables  x.  Because  confidence 
Interval  estimates  and  hypothesis  tests  about  ^ and  P(jc)  are  sensitive  to 


I 


the  variability  in  the  estimated  parameter  vector  attention  has  been 


Inferences  can  be  made  by  the  prediction  of  critical  envelopes  and  by 
testing  relevant  hypotheses.  A critical  envelope  may  be  defined  as  the 
set  of  all  combinations  of  x for  which  the  predicted  probability  of  Injury 
is  less  than  some  given  amount.  However,  because  variability  in  the  pre- 
dicted probability  causes  error  variability  in  the  predicted  critical 
envelope,  this  must  be  taken  into  account  in  the  prediction  process. 

In  addition  to  the  use  of  critical  envelopes,  hypothesis  tests  can  be 
conducted  to  draw  inferences  of  P(x)  in  relation  to  some  small  probability 
of  injury,  P^.  The  statistical  structure  of  hypothesis  testing  can  be  used 
to  control  the  probability  of  making  a wrong  decision  concerning  the  relation- 
ship of  P(x)  to  P^. 

Pertinent  statistical  Inference  procedures  have  been  discussed  in  a 
recent  technical  report  [11].  In  addition,  a Bayesian  approach  to  these 
inferences  has  also  been  explored  [12].  A Bayesian  approach  incorporate."! 
prior  information  with  sample  information  by  the  use  of  a prior  probability 
distribution.  This  distribution  can  be  used  to  quantify  an  experimenter's 
prior  beliefs  and  expert  knowledge  about  an  experiment. 


C.  EVALUATION  OF  MOTION  SICKNESS  PREDICTION  MODELS 

A third  area  of  research  has  been  concerned  with  statistical  problems 
related  to  the  investigation  of  ship  motion  effects  on  humans.  One  goal  of 
this  investigation  is  to  predict  motion  sickness  incidence  (MSI)  associated 
with  complex  motion  characteristic  of  Navy  ships  at  sea.  Other  researchers 
have  defined  MSI  as  the  probability  of  emesis,  and  have  developed  two  MSI 
predictive  models,  but  only  for  sinusoidal  motion. 

Both  of  these  models  appear  to  yield  relatively  accurate  MSI  predictions 
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in  those  cases  where  motion  Is  defined  by  a single  sinusoid.  However, 
neither  model  can  be  used  directly  In  predicting  MSI  for  ship  motion, 
because  that  motion  tends  to  be  broadband.  Desmatlcs  has  examined  three 
approaches  for  obtaining  such  predictions  In  light  of  the  observed  MSI  data 
from  an  experiment  In  which  subjects  were  exposed  for  two  hours  to  motion 
produced  by  the  sum  of  two  sinusoids.  These  three  approaches  may  be  referred 
to  as  (1)  the  Independent  effects  approach,  (2)  the  weighting  approach, 
and  (3)  the  least  squares  weighting  approach. 

The  Independent  effects  approach  assumes  that  when  motion  results  from 

a combination  of  sinusoids,  there  Is  no  Interaction  (l.e.,  synergism)  between 

them  relative  to  MSI.  The  second  approach,  analogous  to  A-welghtlng  In 

acoustics.  Involves  weighting  all  frequencies  back  to  a single  frequency  f^. 

The  frequency  f Is  defined  as  the  "most  critical  frequency",  that  Is,  the 
0 

frequency  which  results  In  the  largest  MSI  prediction. 

The  third  approach,  like  the  previous  one,  also  weights  all  frequencies 
back  to  a single  frequency  f^.  However,  whereas  the  previous  approach 
specifies  the  weights  a priori  based  on  an  existing  sinusoid  model,  the 
empirical  least  squares  weighting  approach  determines  them  a posteriori  from 
experimental  MSI  data  resulting  from  motion  corresponding  to  complex  wave- 
forms. Thus,  this  approach,  unlike  the  other  two,  requires  access  to  MSI 
data  from  motion  that  Is  other  than  sinusoidal. 

Under  the  hypothesis  that  a given  approach  results  In  the  prediction  of 
the  true  MSI  value  (denoted  by  MSIq)  for  a duration  of  two  hours  under  specific 
motion  conditions, 

X - nCMSIg) 


n(MSI  ) (1  - MSr  ) 
0 0 
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would  be  approximately  distributed  as  a standard  normal  random  variable, 
where 

n - the  number  of  subjects  observed  under  this  motion  condition 

and 

X ■>  the  number  of  subjects  who  were  overcome  by  emesis  within 
two  hours. 

Based  on  the  observed  values  of  z calculated  from  the  experimental 
data,  a test  of  the  hypothesis  may  be  made.  However,  when  a number  of 
tests  are  made,  the  problem  of  multi-test  bias  arises.  Multi-test  bias 
refers  to  the  phenomenon  that  as  the  number  of  tests  made  Increases,  the 
chances  of  making  at  least  one  type  I error  also  Increases.  To  avoid  this 
sort  of  contamination.  It  Is  desirable  to  use  a single  overall  test. 

Such  a test  may  be  based  on  the  statistic 
H 

X - -2  Z ln(pj 
1-1  ^ 

where  p^  denotes  the  probability  of  the  observed  outcome  of  the  1— 

Individual  test  In  a total  of  N tests.  Under  the  hypothesis,  the  statistic 
X has  an  approximate  Chi-square  distribution  with  2H  degrees  of  freedom. 

The  results  of  this  test  shed  serious  doubt  on  the  usefulness  of  both  the 
Independent  effects  approach  and  the  weighting  approach  as  methods  of  predict- 
ing MSI.  Tests  of  the  least  squares  weighting  approach  are  Inappropriate 
and  therefore  were  not  performed.  For  that  approach,  extremely  good  agreement 
between  observed  and  predicted  MSI  values  Is  to  be  expected  because  two 
weights  were  estimated  based  on  a total  of  only  three  observations.  Further- 
more, predicted  MSI  values  would  be  compared  with  the  observed  MSI  values 
from  which  the  weights  had  been  estimated.  A valid  test  of  this  approach 
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would  require  a larger  set  of  observed  data. 

Although  only  a small  amount  of  empirical  MSI  data  exists  for  motion 
other  than  that  produced  by  a single  sinusoid,  that  data  provides  a strong 
Indication  that  neither  the  Independent  effects  approach  nor  the  weighting 
approach  produces  accurate  MSI  prediction.  (Judgment  about  the  least 
squares  weighting  approach  must  be  reserved  until  enough  empirical  data 
exists  to  provide  an  adequate  test.)  Details  of  the  statistical  analysis 
of  the  MSI  data  are  given  In  a Desmatlcs  technical  report  [6]. 
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